TeV neutrinos from core collapse supernovae and hypernovae 
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A fraction of core collapse supernovae of type Ib/c are associated with Gamma-ray bursts, which 
are thought to produce highly relativistic jets. Recently, it has been hypothesized that a larger 
fraction of core collapse supernovae produce slower jets, which may contribute to the disruption and 
ejection of the supernova envelope, and explain the unusually energetic hypernovae. We explore the 
TeV neutrino signatures expected from such slower jets, and calculate the expected detection rates 
with upcoming Gigaton Cherenkov experiments. We conclude that individual jetted SNe may be 
detectable from nearby galaxies. 

PACS numbers: 97.60.Bw, 96.40.Tv, 98.70.Sa 



Core collapse supernovae (SNe) are known sources of 
1-10 MeV neutrinos, as detected from SN 1987A Q. Pre- 
viously, 100 TeV neutrinos were predicted from protons 
accelerated in internal shocks in relativisticiets associ- 
ated with gamma-ray burst (GRB) event Q, indepen- 
dently of any SN connection. Long duration (> several 
seconds) GRBs are now known, in several cases, to be 
associated with SN type Ib/c events triggered by a stel- 
lar core collapse which is contemporaneous with the GRB 
event Q . The relativistic jets from such core collapses are 
expected to produce a ~ TeV neutrino precursor burst, 
while the jet is making its way out of the collapsing stel- 
lar progenitor 0, and this may be present even if the 
jet does not manage to burrow through the stellar enve- 
lope, i.e. choked bursts which are v bright, but 7 dark. 
These previous TeV neutrino predictions assumed jets, 
with bulk Lorentz factors Tf, > 100. The frequency of 
highly relativistic jets in SN like events is at most com- 
parable to the ratio of GRB to the average core-collapse 
SN events, < 10~ 3 . 

The deposition of bulk kinetic energy from a jet into 
the stellar envelope may be a powerful contributor to the 
triggering of a core collapse SN or SN Ib/c [5| and/or to 
the disruption of the envelope Q. Only a small fraction 
of core-collapse SNe appear to be associated with GRBs 
• However, a significantly larger fraction of SNe may be 
endowed with slower jets, of T& ~ few 0, SIS El; which 
may also give rise to "orphan" radio afterglows not asso- 
ciated with detected 7-ray emission 11]. In this Letter, 
we point out that SNe endowed with slow jets will be 
detectable sources of multi-GeV to TeV neutrinos. Since 
the occurrence rate per galaxy of SNe is much larger than 
for GRBs, the chances of one occurring nearby is signif- 
icantly enhanced. Hyperenergetic SNe (hypernovae) will 
be fewer in number (~ 10% of type Ib/c rate [l2|). but 
will have a higher neutrino flux. We calculate the neu- 
trino fluence from such slow-jet SNe within the nearest 20 



Mpc, and the corresponding neutrino detectability with 
upcoming kilometer scale Cherenkov detectors in Antarc- 
tica and in the Mediterranian [l^ . 

Jet dynamics. — We take a jet, with r& = 10 r&o.5) 
inside the SN progenitor star. The corresponding jet 
opening angle is 0j ~ 1/Lf,. With a variability time t v = 
0.1i„._i s at the jet's base, internal shocks occur in the 



10.8r2 



r 2 t 

1 fc,0.5 , "' J ! - 1 



jet at a collision radius rj = 2T 2 ct v w 10 
cm, which is inside the progenitor star (typically a Wolf- 
Rayet star of ~ 10 11 cm radius for a SN lb, or larger 
for other type II SN). Typical jet energy, inferred from 
GRBs, is Ej ~ 10 515 i?5i.5 ergs while inside the progen- 
itor star. By analogy with the GRB case, we assume 
that a fraction e e ~ 0.1e ej _i of Ej is converted into rel- 
ativistic electron kinetic energy in the internal shocks, 
which then synchrotron radiate in the presence of mag- 
netic fields representing a fraction Sf, ~ O.le^-i of Ej. 

The volume number density of leptons and baryons 
present in the jet is n' e = n' p = Ej/(2wrjm p c tj) ~ 
10 20 - 5 E 5 i, 5 T^ q 5 t~lt~ 2 _ 1 cm -3 in the comoving frame of 
the jet. Here tj ~ lOi^i s is the typical jet duration. 
Electrons and protons are expected to be accelerated to 
high energy in the internal shocks. Typically all elec- 
trons cool down by synchrotron radiation within the dy- 
namic time scale t' dyn — feu /IV However, due to the large 
Thomson optical depth (r Th ~ W^E^^T^tjlt^), 
synchrotron photons in the jet thermalize to a black body 
temperature of E\ = [15{hc) 3 E J e e /{2TT A r 2 ct J )] 1 ^ « 



4.3£ , s 1 ( 4 5 eef-i r M.5^i V ' i ^-r keV - Also Photons from the 
shocked stellar plasma do not diffuse into the jet as a 
result of the high rxh- The volume number density of 
thermalized photons produced by the shocks is n'^ = 

Ej£ e /(2Trr]cE'tj) 
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The magnetic field strength in the jet is given by B' = 
[AEjS^ctj)}^ » IO^^H-^mV^Vi 
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in the comoving frame. The corresponding shock ac- 
celeration time for protons is t' acc — AE' p /eB' w 

10- u (^/GeV)A 1 ^/ a e 6 -^^ i0 . B ^/ a t„ j _ li s. ^ The 
maximum proton energy is limited by requiring this time 
not to exceed t' dyn or other possible proton cooling pro- 
cess time scales. 

Proton cooling time scales. — The Bethe- Heftier (BH) 
interaction pj —> pe^ has a logarithmically rising cross 
section ctbh = ar2{(28/9)ln[2 J B^/(m p rn e c 4 )]-f06/9}. 
The e^- pairs are produced at rest in the center of mass 
(cm.) frame of the pry collision. The energy lost by the 
proton in each interaction is given by AE' p = 2m e c 2 ^' c m , 
where 7 £. m . = (E' p + E', ( )/{m 2 p c 4 + 2E' p E! y ) 1 / 2 is the 
Lorentz boost factor of the cm. in the comoving frame. 
The energy loss rate of the proton is given by dE' p /dt' BH = 
niyOJBnAEp. The corresponding proton cooling time is 
*bh = E' p /{dE p /dt') = £; p/( 2n 7 CCT BHm e c 2 7 ^ m .) in the 
comoving frame. 

Protons, in the presence of the same magnetic field 
which is responsible for the electron synchrotron losses, 
also suffer synchrotron losses, but over a longer time 
scale given by t' syn — 6Trm p c 3 / (aTh[3 2 mlE p B 2 ) w 
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Protons can also transfer energy to photons by in- 
verse Compton (IC) scattering, because of the sig- 
nificantly high photon density (ni) in the shocks. 
The IC cooling time in the Thomson limit is 
given by t' 1CTh = $m*<? / {^^mlE^E^) « 

3.8(^/GeV)- i ^5^iil1 i0 . 6 t i , 1 ^ i _ 1 s. In the 
Klein-Nishina (KN) limit, the IC cooling time 
is given by ^ CKN = 3E p E! y /(Ao- Th m 2 e c 5 n' 7 ) « 

10- 10 - B (^/(feV)4^6 /a ^i / ll6,0.5*i / i a *«,-l S - The 

Thomson and the KN limits apply for E' p much less or 

greater than m 2 p c 4 /E> 7 « l^E^e^^o^tl^ 
GeV, respectively, in the comoving frame. 

We have plotted the proton acceleration time and the 
different cooling times in the comoving frame as functions 
of the proton energy in Fig. Note that the BH cool- 
ing time (long dashed line) and the synchrotron cooling 
time (short dashed line) are first longer and then shorter 
than the maximum proton acceleration time (thick solid 
line). The corresponding maximum proton energy can 
be roughly estimated, by equating the t' syn to t'^, as 

E' m » l^V^V^ - 6 ^ 4 *^ GeV - 

Pion and muon production and cooling. — High-energy 

protons interact with thermal photons in the shocks 

to produce neutrinos dominantly through j>7 — * A — > 



channel. 
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threshold proton energy at A is E' th = 0.3/ E' 7 pa 
10 4 ' 8 S^" 1 ^ 4 £~i / ' 1 4 r(, i o.5t^ / ' 1 4 t^ 1 GeV and the correspond- 
ing optical depth is t p1 ps ^^E^sl^T^^^^ , 
both evaluated in the comoving frame. However, low en- 
ergy protons may also produce A by interacting with 
photons from high-energy tail of the black body photon 
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FIG. 1: Proton cooling time scales for different processes, in 
the comoving frame, as functions of energy. Also shown is the 
proton shock acceleration time. Burst parameters used are 
E sn = 10 51S erg, tj = fO s, t v = KT 1 s, 17 = 3, e e = e b = 0.1. 



distribution until r p7 falls below unity or below the opti- 
cal depth for other processes such as proton-proton (pp) 
interaction. Similarly, protons with energy above E' p th 
interact with low energy photons to produce A. pp in- 
teractions are effective at lower energy compared to pry 
interactions. The total pion multiplicity from pp inter- 
actions is of the order unity in the energy range we are 
considering here [l4| . For simplicity, we assume that each 
pp interaction produces single pion, the same as pj in- 
teractions, down to the proton energy for which the cor- 
responding neutrino energy is above detection threshold 
(E v ,th) ° n Earth. The available maximum proton energy 
to produce A, through pj interactions, is E' p max « 10 6 ' 3 
GeV. 

High-energy pions and muons produced in pj interac- 
tions do not all decay to neutrinos, as synchrotron ra- 
diation and IC scattering reduce the energy. The syn- 
chrotron break energies (below which pions and muons 
decay before undergoing significant cooling) are found 
by equating their synchrotron cooling times (t' n ^.gyn) 
to their decay times (7^ M idec) in the comoving frame, 



which are E'_ 
and E' „ 
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GeV 

„...,. W^E^e^Tl^X-i GeV, respec- 
tively. The IC cooling time (t' IC ) is comparable to the 
synchrotron cooling time and the two can be combined 
into a single loss rate i c ~ ol ~ + t^ 1 . For IC losses 
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t' lc while in the KN limit 



in the Thomson limit t' 
Kyn ^ *ic> ano - the neutrino fluence in the presence of 
these energy losses is suppressed, compared to the flu- 
ence obtained assuming pions and muons decay without 
any energy loss, by a factor 

O.E-K ,fi) ^7r,/^;cool/^7r./i;dccay ^ ^7T,/^;cool/^7r,/^;syn 

(-^7r,/i/-^7r,/i;sb) 7 E-k^ !> E-x ^-sh. (1) 

Thus muons and their decay neutrinos are severely sup- 
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pressed due to their high energy electromagnetic losses. 
We calculate the neutrino fluence from pion decay next. 

Neutrino flux. — The observed isotropic equivalent 
proton fluence from the SN jet, denoting the quan- 
tities in the observer's frame with subscript "ob", 
at a luminosity distance Dl, is T v , \> = E^T\{\ + 
z)/(2nD 2 L E 2 ^ oh ln[E p>max /Ep !min \). Here we have as- 
sumed a shock accelerated proton energy distribution 
oc Ep 2 . The energy and time in the observers frame 
and in the local rest frame (at red shift z) are related 
by E p = £"p,ob(l + z) and tj — ij j0 b/(l + z). The corre- 
sponding fluence from proton interactions is 
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per SN burst, assuming all shock accelerated protons 
convert to pions. The factor of 1/8 is due to the fact 
that in the absence of energy losses, the neutrinos pro- 
duced by pion decay carry 1/8 of the energy lost by pro- 
tons to pion production (since charged pions and neu- 
tral pions are produced with roughly equal probability, 
and since the muon neutrino produced in a pion de- 
cay carries roughly 1/4 of the pion energy). The neu- 
trino energy range is then -E„ i0 b,min - E v 

,ob,max 

- 0.05T b E'/(l + z 
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2.5 105-5 



10 5 - 5 /(l + z) GeV in 
the observer's frame on Earth. Here we have assumed 
Eujh. = 10 2 ' 5 -Ev,2.5 GeV for a typical ice Cherenkov 
detector such as IceCube. For a SN at ~ 20 Mpc 
(Dl ~ 10 25 - 8 cm, z ~ 0), e.g. in the Virgo cluster, the 
neutrino fluence would be 
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GeV^cm- 2 ; 10 2 ' 5 < £„, ob /GeV < 10 5 5 .(3) 

To calculate the diffuse neutrino flux, we sum over flu- 
ences from all slow-jet SNe distributed over cosmological 
distances in Hubble time. The SNe rate follows closely 
the star formation rate (SFR) which can be modeled 
as p*(z) = 0.32/i 70 exp(3.4z)/[exp(3.8z) + 45] M yr" 1 
Mpc -3 per unit comoving volume [15j. Here Hq — 70ft.7o 
km s" 1 Mpc" 1 is the Hubble constant. For a Friedmann- 
Robertson- Walker universe, the comoving volume ele- 
ment is dV/dz — 4ttD l c/(1 + z)\dt/dz\ and the rela- 
tion between z and the cosmic time t is (dt/dz)^ 1 = 
-H (l + z)y/(l + tt m z)(l + z) 2 -ft A (2z + z 2 ). Here we 
have used the standard ACDM cosmology with O m = 0.3 
and f^A = 0.7. 

The number of SNe per unit star forming mass (/ sn ) de- 
pends on the initial mass function and the SN threshold 
for stellar mass (M ~ 8 M Q ). A Salpeter model <j)(M) oc 
M~ Q with different power-law indices can generate dif- 
ferent values for / sn , e.g. « 0.0122 M" 1 for a = 2.35 [3 



and w 8 x 10~ 3 M" 1 for a = 1.35 0. We adopt the 
model in Ref. 0] which corresponds to the local type 
II SNe rate n sn (z = 0) = f sn p*(z = 0) 
yr -1 Mpc -3 agreeing with data [T^ . 

The distribution of SNe per unit cosmic time t and 
solid angle O covered on the sky can be written as 
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We assume a fraction £ sn < 1 of all SNe involve jets and 
a fraction l/2r 2 of all such jets are pointing towards us. 
The observed diffuse SNe neutrino flux, using Eq. J5J), is 
then 
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We have plotted the diffuse flux from all cosmologi- 
cal slow-jet SNe in Fig. Elby numerically integrating Eq. 

assuming the maximal fraction £ sn = 1. Also shown 
are the cosmic ray bounds (WB limits) on the diffuse 
neutrino flux 0], and the atmospheric v u , v, L flux from 
pion and kaon decays (conventional flux, 18]), compati- 
ble with AMANDA data [lflj . 
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Figure |3 demonstrates that the diffuse flux is unlikely to 
be detectable by km-scale telescopes. However, individ- 
ual nearby SNe may be detectable, as we show below. 

Neutrino event rate. — The effective area (-A c ff) of a 
Cherenkov detector depends on the zenith angle and the 
energy of the muon that is created by a charged cur- 
rent neutrino-nucleon (vN) interaction in the vicinity of 
the detector. IceCube can trigger on muons of energy 
(which is ~ 80% of the neutrino energy ) ab ove a few hun- 
dred GeV created within A c $ ~ km 2 [13|. We assumed 
E v ,tb. ~ 300 GeV, as discussed before. The angular sensi- 
tivity of the IceCube detector is 1° for muon tracks with 
< 140° zenith angles below 1 TeV. The pointing resolu- 
tion gets better at higher energy. 

The likeliest prospect for detection is from individ- 
ual SN in nearby starburst galaxies, such as M82 and 
NGC253 (Dl = 3.2 and 2.5 Mpc in the Northern and 
Southern sky, respectively), over a negligible background, 
using temporal and positional coincidences with optical 
detections. The SN rate in these is ~ 0.1 yr" 1 [2(J, much 
larger than in our own galaxy or in the Magellanic clouds. 
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FIG. 2: Diffuse muon neutrino flux on Earth from all (£ sn = 1, 
Ej = 10 51 ' erg) core collapse SNe (heavy solid curve). The 
dashed curves are the cosmic ray (WB) limits, and the light 
full curve is the atmospheric v^. and flux. This shows that 
the diffuse SNe flux is unlikely to be detectable with km 2 de- 
tectors. However, individual nearby SNe may be detectable, 
as discussed below. 

To calculate the number of events in a km 2 detector 
such as IceCube, we use the fluence of Eq.© from each 
SN in the 300 GeV-300 TeV energy range, and the full 
detection probability depending on the source's angular 
position. This probability depends on the Earth's shad- 
owing effect and the energy dependent uN cross section 
inside A e s |2lJ. The number of events from M82 and 

NGC253 is 281 E^[ 2 5 and 482 respectively per in- 

dividual SN with its jet pointing towards us. The timing 
uncertainty in the optical detection of SN is ~ 1 day. The 
corresponding atmospheric background events within 1° 
angular resolution is 0.07 for both M82 and NGC253 in 
the same energy range. At the quoted rate, a SN from 
one of these galaxies would be expected within five years. 
Other nearby spirals (M31, M74, M51, M101, etc, and 
the Virgo cluster) will also contribute. At the standard 
rate of 1 SNu= 10 _2 yr _1 per 10 10 blue solar luminos- 
ity for average galaxies, the roughly 4000 galaxies known 
within 20 Mpc would suggest a rate > 1 SN yr _1 . Ice- 
Cube may detect roughly 1/5 of them at a level of > 1.5 
Up events/SN. 

Because of oscillations, neutrinos of all three flavors 
(up, u e , u T ) should reach us in equal proportion. How- 
ever, only u^'s are emitted from the sources under con- 
sideration and the total number of neutrino events (in- 
cluding v e and v T ) will remain the same as the events 
we have calculated. The lack of good directional sen- 
sitivity for the v e and v T events may prevent obtaining 
their positional coincidence with the SN. (Also, the A e s 
for v e and u T detection are somewhat smaller than u^.) 
However the timing coincidence of u e and v T events with 
Up events may still be useful to verify the neutrino oscil- 
lations at these energies, and test their common origin. 

Discussion. — While ultra relativistic jets are thought 



to be responsible for the long-duration GRB associated 
with a small fraction of massive core collapse SNe, mildly 
relativistic jets may occur in a much larger (< 1) fraction 
of core collapse SNe. Such slow jets may contribute to the 
bounce needed for ejecting the stellar envelope, resulting 
in the observed SN optical display. Late time (~ yr) 
radio emission may be a signature of slow jets [llj . Such 
jets may also be needed to explain the apparent unusually 
large energies of hypernovae and the anisotropics inferred 
from optical polarization measurements of core collapse 
SNe. Here we propose, as an independent test of this 
model, the observation of TeV neutrinos, which may be 
detected in the near future with IceCube and other km 2 
neutrino Cherenkov detectors. 
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